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At hole concentrations below x=0A, Bai-^K^BiOa 
is non-metallic. At x — 0, pure BaBiOs is a Peierls 
insulator. Very dilute holes create bipolaronic point 
defects in the Peierls order parameter. Here we find 
that the Rice- Sneddon version of Peierls theory pre- 
dicts that more concentrated holes should form stack- 
ing faults (two-dimensional topological defects, called 
slices) in the Peierls order parameter. However, the 
long-range Coulomb interaction, left out of the Rice- 
Sneddon model, destabilizes slices in favor of point 
bipolarons at low concentrations, leaving a window 
near 30% doping where the sliced state is marginally 
stable. 
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I. INTRODUCTION 

When a ID crystal is driven incommensurate, 
for example, by doping to alter the Fermi wavevec- 
tor, it is well-known that the modulation of the 
commensurate crystalline order tends not to be 
uniform but to accumulate at kinks in the crys- 
talline order parameter 10 . A prototype for this is 
the Su-Schrieffer-Heeger model for polyacety- 
lene, (CH) X . In a simple chemical view, the car- 
bon chain has alternating single and double bonds. 
Crystallography confirms this lattice dimerization. 
The scalar order parameter (amplitude of the stag- 
gered charge disproportionation) has two possible 
values, ±p. Zero-dimensional domain walls (topo- 
logical defects) separate regions of the ID chain 
with positive and negative order parameters. Su et 
al. showed that when carriers are doped into poly- 
acetylene, they create new domain-wall defects and 
localize into mid-gap electronic soliton states on 
these defects. Analogous effects in two dimensions 
H occur in layered materials near a Mott insula- 
tor phase, such as cuprates, nickelates, and layered 
manganites. The one-dimensional domain wall in 
the 2D system is generally called a "stripe." The 
3D analog, which occurs in some cubic manganites 
B , is a planar domain wall || , sometimes called a 
"sheet" or "lamella." We prefer the term "slice," 
which, like the word "stripe" can be used as verb or 
noun. Here we discuss a 3D model system in which 
bonding and strain energies favor slices, but long- 
range Coulomb interactions prefer point defects. 
We analyze the system numerically and show that 



at doping levels near 30%, Bai-^K^BiOa is near 
the crossover where slices may become more favor- 
able than point defects for lowering the energy of 
the doped-in holes. 

Charge-density-wave (CDW) and Peierls sys- 
tems have scalar order parameters pq = 

N^ 1 J2e Pi ex p(*Q • BaBiOs is a simple 3D ex- 
ample. With nominal valence Bi 4+ , the Bi 6s-band 
is half-filled. The material is non-metallic, with 
a 2eV optical gap 0. Crystallography M shows 
a doubled unit cell; pi alternates with wavevec- 
tor Q = (7r,7r,7r). The simple cubic sublattice 
of nominal Bi 4+ ions self-organizes into a bipar- 
tite (rocksalt-type) charge-ordered array of nomi- 
nal Bi 3+ (p = 4 — pq, called "A" sites) and nomi- 
nal Bi 5 + (p = 4 + p Q , called "B" sites) ions. The 
actual value of the order parameter pq has magni- 
tude < 1, and takes two degenerate values, ±|pq|- 
A simple way to think of this is as a 3D version of a 
Peierls instability. In a bi-partite lattice with only 
nearest-neighbor hopping there is accidental Fermi 
surface nesting: e(k) = — e(k + (ir, ir, ir)), and both 
are zero at the Fermi energy. This guarantees that 
the crystal can reduce its electronic band energy 
via the electron-phonon interaction by dimerizing. 
A Hamiltonian which contains this effect was in- 
troduced by Rice and Sneddon 




FIG. 1. Fragment of BaBi03, showing Bi atoms 
only. Filled circles are Bi 3+ (A) ions, and empty cir- 
cles are Bi 5+ (B) ions. In (a) the normal ABABAB... 
sequence of (111) layers is shown; in (b), a hole-type 
stacking fault ABBABA..., or slice, is shown. 

Regions of charge-ordered BaBiC>3 with pq > 
are separated from regions with pq < by stacking 
faults. The simplest stacking fault lies in a (111) 
plane. In perfectly ordered BaBiC>3, (111) planes 
are alternating A and B type (Fig. |l|a). A stack- 
ing fault with no nuclear disorder has either two 
adjacent A layers (local charge excess -1 per site 
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on the plane, or electron-doped) or two adjacent B 
layers (local charge excess 1, or hole-doped) as in 
Fig. 0b. 

Here we point out that the Rice-Sneddon Hamil- 
tonian, given below, predicts that holes or elec- 
trons, when introduced, will self-organize into 
slices. However, the long-range Coulomb inter- 
action, neglected in the Rice-Sneddon model, will 
destabilize slices except possibly at the most favor- 
able doping level. For light doping, the preferred 
structure for holes is to self-organize into point 
defects, small bipolarons, as previously discussed 




II. RICE-SNEDDON MODEL 

The Rice-Sneddon model || is simple, well- 
studied §|^,H|Lf], and quite successful ||@. 
In the perovskite crystal structure of BaBiOs, each 
Bi atom is surrounded by six oxygens, and each 
oxygen is shared by two Bi atoms. At high tem- 
peratures (of order 1000K) the crystal is nearly 
cubic perovskite, but at lower T there are rota- 
tions and distortions of the BiOg octahedra, which 
enlarge the unit cell. We believe that the most 
important effect is the "breathing" displacements 
u as ±0.12A j7|, p"7||lS| ] of oxygens away from or 
toward alternate Bi atoms. These provide a natu- 
ral electron-phonon mechanism to enhance the Bi 
6s charge density on A sites where oxygens breath 
outward, and reduce the Bi 6s charge density on 
B sites where the oxygens breathe in. A micro- 
scopic Hamiltonian containing the minimal neces- 
sary electron-phonon interaction was given by Rice 
and Sneddon pL 

w = - 1 c \ c v - 9j2 e (t) c W 

<e,e'> e 

+ 1 -KY,u(i+^) 2 . (1) 

(.a 

The first term is nearest-neighbor hopping of Bi 6s 
electrons with hopping integral t as 0.35 eV. The 
index of summation £ implicitly includes a spin as 
well as site quantum number. The filling is 1 — x 
electrons per site. The variable u(£ + a/2) (with 
a = x,y, z) is the displacement along a Bi-O-Bi 
bond in the a direction of the oxygen located at 
position (£ + a/2)a. The variable e(£) is the local 
dilation or "breathing" amplitude of the 6 oxygens 
which surround the Bi ion at site £, 

e(£) = [u{£, x + x/2) - u{$- x/2)] + [x -> y] 
+ [*-«]. (2) 

The Einstein restoring force K as 19 eV/A 2 is 
fitted to the measured 70 meV frequency of the 
Raman-active Peierls breathing mode M] . 



At half-filling, this model opens a Peierls gap 
at the Fermi level; the electron-phonon interaction 
parameter g as 1.39 eV/A is fitted to the measured 
H as2 eV gap. This is an ordinary size electron- 
phonon coupling. The change in Coulomb field of a 
charge -2e oxygen ion gives g as 1.2 eV/A. The re- 
sulting dimensionless coupling constant T = g 2 / Kt 
is «0.30, intermediate between the weak (r < 0.2) 
and strong (r > 0.4) coupling regimes. In this mid- 
dle regime, neither hopping nor electron-phonon 
energy dominates pl|. The ground-state of un- 
doped BaBiOs is as close to a bipolaronic crystal 
(large T, \pq\ as 1) as to the conventional Pcicrls- 
CDW (small T). We calculated §J| the order pa- 
rameter pq at r = 0.3 to be 0.82. The correspond- 
ing oxygen displacement uq — 2gpQ/K — 0.12A 
agrees with the diffraction and EXAFS measure- 
ments 0,0^8|, showing that the model is inter- 
nally consistent. 

We previously found that excitations across the 
Peierls gap form self-trapped excitons |D|. We 
also reported that holes inserted into BaBiOa self- 
trap and form bipolarons jl^] since the coupling 
strength exceeds r c = 0.17. These are doubly 
charged point defects, corresponding to local de- 
pressions of the order parameter where the oxy- 
gen distortion e(£) — > for t = 0. For non-zero 
hopping t the bipolaron spreads out and evolves 
continuously from a small bipolaron (r ^> T c ) to 
large CDW-like bipolaron as T — > r c . The stability 
of bipolaron defects provides a simple explanation 
why dilutely doped BaBiOa remains insulating and 
diamagnetic. 

The "disproportiona- 
tion reaction" 2Bi 4+ — >Bi 3+ +Bi 5+ has been much 
discussed in the literature on BaBiOs. Using the 
definition 2U = £(Bi 3 +) + £(Bi 5 +) - 2£(Bi 4+ ), 
one can say that the effective Hubbard U parame- 
ter is negative. Two factors may contribute to U, 
electron-phonon effects (expected to be attractive) 
and Coulomb repulsion (U = U cp + U e \). If one 
wants to assign the mechanism for disproportiona- 
tion completely to Coulomb interactions, then the 
Hubbard U calculated with all atoms held station- 
ary in cubic perovskite positions (defined as U e \) 
should be negative. Vielsack and Weber |l(| did 
careful calculations of U c \, finding no evidence for 
negative values, but instead a small positive value 
U c \ as 0.6 ± 0.4eV. Therefore we shall temporarily 
ignore the on-site Coulomb repulsion U c \. Appar- 
ently the Bi 4+ ion must reorganize its environment 
in order to stabilize the charge disproportionation. 

III. HOLE DOPING 

What happens at finite doping concentrations 
x? Assuming sufficient electron-phonon coupling 
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to destroy the undistorted metal, there are two 
possibilities, (1) bipolarons, and (2) slices. 

(1) Numerical studies by Yu et al. |To| , and con- 
firmed by us, show that randomly located point 
bipolarons are at least metastable. Depending on 
whether bipolarons attract or repel, the system 
could then either phase-separate into undoped and 
doped regions, or form spatially separated bipo- 
larons. When bipolarons are small, the energy 
of an array of bipolarons is described by a pair- 
wise additive potential V(T,r), containing the re- 
pulsive long-range Coulomb interaction Vcoui (ne- 
glected for the time being) and the interaction Vt, 
induced by hopping. Vt decays exponentially with 
the distance r between two bipolarons (like bipo- 
laron wave functions.) Since bipolarons can only 
sit on former A sites, the nearest neighbor interac- 
tion Vb = Vt(r, V2a), which is also the strongest 
interaction, is between bipolarons separated by 
a(110).) We computed Vq numerically by optimiz- 
ing the oxygen-positions self-consistently for given 
bipolaron positions. Fig. || shows that large CDW- 
like bipolarons attract each other, whereas small 
bipolarons repel, with Vt — > in the atomic limit 
t — > 0. Perturbation theory around t = gives 
a bipolaron repulsion. At the physically relevant 
value T = 0.3, bipolarons repel according to Fig. 
|^, but multi-bipolaron interactions become impor- 
tant with decreasing T. Therefore, in the interme- 
diate coupling regime we must rely on exact nu- 
merical diagonalization of Eq. (nl) |10| . 



average hole charge on phase-slip B-sites is thus 
phole ~ + 5 ■ The actual values of hole charge found 
for r = 0.3 and x = 1/4 are shown in Fig. ||. 




FIG. 2. Interaction potential Vo for two bipolarons 
sitting on (0,0,0) and (o/2)(l, 1, 0). Large bipolarons 
attract whereas small bipolarons repel. 

(2) In contrast to bipolaron defects where the 
order parameter never changes sign, holes could 
form topological defects. Phase-slips are planar 
defects with sign-changes of the charge order pa- 
rameter pq and of the breathing order parameter 
e = (— 1) e;. Consider a "BB" stacking fault in 
the Ill-direction (Fig. lb). In the atomic (t — > 0) 
limit, each B-site on a phase-slip plane has three 
displaced and three undisplaced oxygen neighbors, 
i.e. ei w ±3mo- Phase-slips accumulate one hole 
for every two atoms on a (111) BB bilayer. The 
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FIG. 3. Normalized breathing order parameter 
e/6uo (upper panel) and hole charge phoic (lower 
panel) for a x = 1/4-doped phase-slip struc- 
ture (AB*B*ABAB*B*)jv at T = 0.3 (two periods 
shown). The order parameter changes sign across each 
phase-slip B*B*. Holes are localized on phase-slip 
planes with slight spreading to neighboring A-planes. 
The hole charge is computed relative to a reference 
system with the same structure, (AB*B* ABAB*B*)jv, 
but with A charges fixed at their undoped value 4+pq, 
B charges fixed at 4-pQ, and B* charges fixed at 4~pq/2 
corresponding to zero doping, with the charge deficit 
Pq spread equally on the two adjacent B* layers. The 
charge fluctuation pQ has the value 0.82 at T = 0.3. 

Our aim is to find the ground-state hole arrange- 
ment, testing the stability of several bipolaronic (1) 
versus the phase-slip (2) solutions numerically. We 
did a series of calculations oni= 1/4-doped sys- 
tems. The bipolaron sytems (1) were (a) maximal 
spacing between bipolarons, obtained when they 
occupy center and corner sites of a tetragonal 16 
atom unit cell (bet structure); (b) a simple cubic 
(sc) arrangement of bipolarons sitting at the corner 
sites of a 8 atom cubic cell; (c) a disordered struc- 
ture with random bipolaron positions; (d) phase- 
separated structures based on unit cells containing 
8 or 16 (111) planes, where we replace one or two 
near neighbor A planes with bipolaronic B planes, 
[(AB) 3 B 2 ]at or [(AB) 6 B 4 ]jv /2 . We also looked at 
unit cells with 8 or 16 111 planes containing phase 
slips (2) [AB (ABB) 2 ] at and [(AB) 2 (ABB) 4 ]a7 2 . 
Finally, we looked at the undistorted metal [C]jv 
where each atom C has a nominal Bi 4 25+ valence. 

Phase separation was strongly disfavored, while 
separated bipolarons and phase slips were all 
metastable. The phase slips weakly repelled, pre- 
ferring the 8 plane solution to the 16 plane solu- 
tion. The order parameter and hole charge density 
of this 8 plane solution are shown in Fig.(^). The 
stability is determined computing the total energy 
-Etot({uj-}) given by Eq.(Q), which is a function of 
the oxygen-positions ur. We start by guessing oxy- 
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gen positions to get close to a local minimum in the 
energy landscape, then vary oxygen positions us- 
ing a gradient minimization routine to find a self- 
consisistent minimum. For smaller periodic struc- 
ture, we used fc-space sampling in the correspond- 
ing Brillouin zones (8000 k-points). For each k, 7i 
is diagonalized exactly. States arc filled with two 
electrons up to the desired doping. For the ran- 
dom bipolaron structure we used large asymmet- 
ric clusters («400 atoms) with periodic boundary 
conditions and k = only. Initial oxygen-positions 
had Peierls order with small random deviations. 
Our calculations on random configurations gener- 
ally reproduce the earlier calculations of Yu et al. 



but we believe that such corrections are not the 
relevant ones. There are also (i) small structural 
distortions (rotations of oxygen octahedra) beyond 
the breathing-mode distortions considered here, 
and (ii) non-adiabatic effects (such as zero point 
motion) associated with the fact that the oxygen 
mass is not infinitely large compared with the elec- 
tron mass. We believe that both of these also have 
little relevance. It is harder to dismiss two other 
effects: (iii) the disorder caused by the dopant 
atoms, and (iv) the long-range Coulomb interac- 
tion, both omitted so far. Of these, the last is 
clearly important, as we now show, and tends to 
destabilize phase-slips relative to distributed bipo- 
larons. 
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FIG. 4. Energy difference E/t of bipolaron struc- 
tures (la-d) relative to the most favorable phase-slip 
structure (2) for x = 1/4-doping. 

Energies of various states at doping a;=l/4 are 
shown in Fig. ^. In the atomic limit t = all bipo- 
laron and phase-slip structures are degenerate. All 
that is required is that no A site should have an A 
first neighbor. Below a critical coupling strength 
T c (x) (ps 0.2 for doping a;=l/4), distorted struc- 
tures become unstable with respect to the undis- 
torted metallic structure. Above r c (a;), we find 
numerically that at x — 1/4 holes strongly prefer 
to order in phase-slips (2). Bipolaronic structures 
(1) are quite similar to each other in energy, and 
behave as expected from the the bipolaron-pair in- 
teraction. In the intermediate coupling range there 
occurs a cross-over from the tetragonal to the lay- 
ered BBB-structure as stable bipolaron configura- 
tion, corresponding to a change in the overall bipo- 
laron interaction from weakly repulsive to weakly 
attractive. For 1/4-doped BaBiOa, the energy gain 
for the phase-slip solution compared to the bipo- 
larons is about 50-65 meV per hole. 

Thus, contrary to previous studies of doping of 
this model |9[jl0|], we find that the stable dop- 
ing state is not bipolarons but phase slips. Does 
this model correspond sufficiently to reality for 
BaBiOs? Iwano and Nasu (l4| use a more com- 
plicated hopping and electron-phonon interaction, 



IV. LONG-RANGE COULOMB EFFECT 

At low doping, there is a large Coulomb 
cost in putting charges onto stacking faults in- 
stead of widely distributed point charges. We 
modeled this as follows. The Madelung en- 
ergy was computed by the Ewald method for 
(Ba 2+ ) 2 (Bi 3 +Bi 5 +)(0 2 -) 6 . The calculation was 
re-done for many large unit cells, with holes added 
on selected Bi ions (Bi 3+ -^Bi 5+ ), and com- 
pensating negative charges distributed uniformly 
throughout space. This is one way to mimic potas- 
sium doping of BaBiO 3 . First consider the sliced 
state. By numerical calculation for uniformly dis- 
tributed phase slips at many values of x between 
1/39 and 1/3, we found a good fit to the for- 
mula £5 = (-25.1 + n/9x)(e 2 /2ae 00 ). The static 
electronic screening ps 5 was measured by 
Tajima et al. [^0|. Eg is the difference of energy 
per hole between the sliced solution and the un- 
doped Peierls insulator. The term 7r/9ir is the an- 
alytic result for idealized uniform sheets of charge 
of vanishing thickness, arranged periodically in a 
compensating charge background. The term -25.1 
corrects for the discreteness of the charges, the 
absence of the self-interaction, and includes the 
Coulomb energy of the holes with the background 
BaBiOa lattice. 

We also need the energy difference per hole Eb 
of optimally spaced bipolarons relative to the un- 
doped Peierls insulator. For doping x — 1/n 3 the 
bipolarons can be placed on a sublattice of face- 
centered cubic (fee) form, with maximum spacing. 
Numerical results for n—2, 3, and 4 fitted well to 
the formula E B = f— 21.2 - 4.585a; 1 / 3 )(e 2 /2ae 00 ). 
The term —4.585a; 1 ' 3 is the Madelung energy of an 
fee lattice of charges 2e in a uniform compensating 
background, and the constant -21.2 accounts for 
the energy of the holes in the background BaBiOs 
lattice. The difference Eg — Eb is the Coulomb 
penalty per hole for forming slices or phase slips. 
It is plotted in Fig o, becoming large at low x with 
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a shallow minimum near x = 0.3. 

A more accurate estimate of the Coulomb 
penalty is shown as an x on Fig and was ob- 
tained by doing Ewald sums using the actual Bi 
site charges (pi = 5 — rii where rij = 2 < ipjipi > 
is the local occupancy of the Bi s states) and 
the actual (AB*B*ABAB*B*)at slice structure and 
bet bipolaron structure. It is seen on the fig- 
ure that at 1/4 doping, the Coulomb penalty is 
twice as high as the Peierls benefit in forming 
slices. These calculations did not account for the 
actual random positions of the compensating neg- 
ative charges where Ba 2+ ions are replaced by 
K + ions. Instead, the compensating charge was 
distributed uniformly. The calculations were re- 
peated for a model where the compensating nega- 
tive charges were equally shared by each Ba atom 
(Ba 2+ -^B&( 2 ~ x ^ + ). This only changed the con- 
stant terms in E$ and Eb, but did not affect the 
difference, Es — Eb plotted here. 



Weber [JlfjJ, the onsite correction 0.07 ± 0.05eV is 
potentially sufficient to re-stabilize the sliced solu- 
tion. Thus we can say that the best place to look 
for sliced structures in BaBi03 is in the range near 
1/4 to 1/3 doping, but we cannot predict whether 
the sliced solution will be destroyed by disorder or 
Coulomb effects. 

Doping 1/3 is the highest at which a simple 
sliced solution ((ABB)n) is possible. For higher 
doping, either the A-planes acquire polaron defects 
or the sliced solution is destroyed and purely bipo- 
laronic states win. 

It is worth mentioning that there have been 
reports |H] (subsequently attributed to electron 
beam heating effects [^2|) of superlattice diffrac- 
tion spots in doped BaBi03, not apparently iden- 
tical to the superstructures predicted here. A fur- 
ther search in the doping region 0.2 < x < 0.35 
would be interesting. 
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FIG. 5. The curve is an analytic approximation to 
the Madelung or Coulomb repulsion of a periodic ar- 
ray of (111) phase-slips (slices) compared with opti- 
mally spaced separated bipolarons at the same av- 
erage density, where doped-in holes are modelled as 
point charges (charge =e) with screening set to 
5. The open circles give the numerically computed 
Peierls short-range stabilization energy (from hopping 
and electron-phonon energy) of sliced solutions relative 
to separated bipolarons. The symbol x is the exact nu- 
merical Madelung energy difference at 1/4 doping be- 
tween the sliced solution and the bet bipolaron solution 
with the numerically computed layer charges pi. 

The Madelung sums discussed above are not the 
complete Coulomb effect. The missing on-site re- 
pulsion has the opposite effect, preferring sliced 
solutions to distributed bipolarons. If we add back 
the local term 



(3) 



and treat it as a first-order perturbation, the sliced 
solution at 1 /4-doping has lower on-site energy by 
0.121C/ei per hole than the bet bipolaron solution. 
Using the value U c \ « 0.6±0.4eV from Vielsack and 
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